Sum rules for spin-Hall conductivity cancelation 
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It has been shown recently that the universal dc spin conductivity of two-dimensional electrons 
with a Rashba spin-orbit interaction is canceled by vertex corrections in a weak scattering regime. 
We prove that the zero bulk spin conductivity is an intrinsic property of the free-electron Hamiltonian 
and scattering is merely a tool to reveal this property in terms of the diagrammatic technique. When 
Zeeman energy is neglected, the zero dc conductivity persists in a magnetic field. Spin conductivity 
increases resonantly at the cyclotron frequency and then decays towards the universal value. 

PACS numbers: 72.25.-b 
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In the framework of the burgeoning field of semicon- 
ductor spintronics ,"tS'ij i an active interest currently exists 
in the mechanism of generating spin fluxes by electric 
fields. Pure spin fluxes that are not accompanied by elec- 
tric currents are expected to reduce dissipative losses re- 
lated to spin injection and spin transport. Therefore, 
the proposals by Murakami et al.'^ and Sinova et alA 
of generating dissipationless tranverse spin currents (a 
spin-Hall effect) by a driving dc electric field acquired 
a lot of attention. For two-dimensional (2D) systems 
with a spin-split (Rashba) energy spectrum the univer- 
sal spin conductivity of Ref. |^ has been put under thor- 
ough scrutiny. It turned out that the dissipationless com- 
ponent of the bulk spin current, which originates from 
virtual inter branch transitions, is canceled by ladder di- 
agrams when a weak impurity scattering is taken into 
account. This cancelation was originally discovered by 
Schwab and Raimondi^ in a somewhat different context, 
rediscovered by Inoue et al.J and recently confirmed in a 
number of papers i^i^iiSiii The cancelation is rather puz- 
zling because it does not follow from any known symme- 
try arguments. Therefore, it is important to uncover its 
origin. 

The Hamiltonian in question is 



H = h^k /2m + a{cr x fe) • z. 



(1) 



and in the absence of magnetic field, B = 0, spin currents 
are defined as 



(2^ 



{\\a,Vj{k)+Vj{k)a,\X) (2) 



with a proper choice of the integration area in the k- 
space. Here v — h^^dH /dk is the velocity opera- 
tor, j = X, y, and z is a unit vector perpendicular to 
the confinement plane. Spin currents J^zx and J^zy are 
driven by electric fields Ey and E^, respectively. The 
corresponding spin conductivities are Y,zxy = Jzx/Ey 
and Yizyx = JzylEx- Because the Hamiltonian H pos- 
sesses a Coov symmetry, applying (Ti, reflections results in 



-'zxy 



-Yjzyx- Hence, the part of the tensor related 



to the two last indices is equivalent to M^, a z compo- 
nent of a pseudovector M orthogonal to the xy plane, 
and a product OzMz is an invariant of the group Coov 



These arguments prove that T^zij is an invariant tensor 
of the spatial group. It also invariant with respect to the 
i-inversion. Under these conditions, one can expect that 
vanishing of S^y is a result of some sum rules, and we 
present these sum rules in what follows. 

The oscillator strength of a free electron state related 
to the "transition to itself manifests itself in the absorp- 
tion as a cyclotron or Drude spectrum, i.e., only in the 
presence of a perturbation like a magnetic field or some 
electron scattering mechanism. The same problem can be 
relevant for spin currents, hence, we impose a magnetic 
field B II z by considering k in Eq. as a kinetic mo- 
mentum k = — iV + eA/?ic, but do not include a Zeeman 
term into the model Hamiltonian H . The exact solution 
of this problem is^^ 



Exn = huJcin + Ac„) (3) 



where 



V2 



y/-f^n + 1/4, A = -A. (4) 



Here 7 = [2{ma'^ /h^)/hLJc]^/^ is a dimensionless spin- 
orbit coupling constant, tuvc — eB/Tic is the cyclotron 
frequency, A = ±1 designates the spectrum branches, n > 
1, and ipn are oscillator wave functions. The eigenstate 
n = exists only for A = 1, with bQ = l and 6j = 0. 

The standard Kubo-Greenwood formula results in the 
following expression for the frequency dependence of the 
real part of spin conductivity 



'^zxyi^) 



^ {Xn\vy\X'n'){X'n'\Jzx\Xn) 
/ . 7~. ~. vi ^2 ■ y^) 
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[uJx'n' - LOXnY 



The asterisk over the summation sign indicates that sum- 
mation should be performed only over the states with 
Loxn = Exn/Ti > m/^j tox'n' < A*/^: M being the Fermi 
energy; below, we accept ^ > 0. Here / — {ch/eBy^"^ 
is the magnetic length. It appeared through the fac- 
tor 1/2ttP accounting for the Landau level degeneracy. 
When deriving Eq. (O we took into account that the ma- 
trices {Xn\vy\X'n') and {Xn\J'zx\X'n') are antisymmetric 
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and symmetric, respectively. For the same reason, the 
sum includes only nondiagonal matrix elements. 

It is seen from Fig. 1 that for A 7^ A' the summation 
over n should be performed over a vast area, while for 
A = A' the sum includes, for each branch, only a single 
term coming from the two states adjacent to the Fermi 
level. It will be shown in what follows that these two 
terms are anomalously large and in the dc regime cancel 
the sum. 

To simplify Eq. one can take advantage of the 
fact that both operators in the numerator can be ex- 
pressed through commutators including the Hamiltonian 
H. First, Vy = {i/h)[H,y]-, where the operator of the 
coordinate y = yo — Pkx, and yo is the center-of-orbit 
operator^ The operators kx and ky are related, in the 
standard way, to the creation and annihilation operators, 
kx — {a^ + a)/ 1^/2 and ky = {a'^ — a)/il^/2. Because 
commutes with a and a"*", it can be disregarded in non- 
diagonal matrix elements of Eq. lO, hence, 

{\n\vy\X'n') = -ifiujxn - LLi\'n'){\n\kx\>^ n') . (6) 

Second, the spin-current operator Jzx equals Jzx = 
— {ifi/2ma)[H, ax]- as has been noticed by Dimitrova«ii 
Therefore 

{\'n'\Jzx\Xn) = {uj\'n' - u!xn){>^'n'\ax\Xn). (7) 

Zma 

After substituting Eqs. © and Q into Eq. (|SJ), we arrive 
at a simple equation for "Szxy = '^zxyi^ — 0) 



because k- — k+ = 2ma/'h^ . This result coincides with 
the universal conductivity by Sinova et al^ (with the 
accuracy to the factor h/2 due to the difference in the 
definition of J^zx)- 

Equation (|10|l indicates that the interband contribu- 
tion to '^zxyi that includes a sum over the whole "bulk" 
of the Fermi sea, in the B ^ Q limit is expressed in 
terms of the Fermi surface parameters, the momenta fc_ 
and This result resembles the recent inference by 
Haldanai^ that the anomalous Hall effect is a Fermi sur- 
face property. 

Two intrabranch contributions to T,zxy come from the 
diagonal in A terms in Eq. © 



zxy 



ieh 
27rma 

Xehj {nx-l)cn^+n\Cn>,-i + X/2 



8Tv^/2mal 



{Xnx\kx\Xnx - l){Xnx - l\ax\Xnx) 

(11) 



Here nx is the lower empty orbit at the A branch, A = ±1, 
hence, these contributions are Fermi surface parameters. 
For a weak magnetic field, Eq. reduces to 



XX _ X ehkx 



Therefore, in the B ^ limit. 



^zxy 



ieh 
2'nma 



^zxy — '^zxy 



^zxy 



^zxy ^' 



(12) 



(13) 



^ (An|fc^|AV)(A'n'|cr^|An). (8) in agreement with the conclusions of Refs.|Hi|lT3 We 



AnA' 



It is seen from Fig. 1 that there are two types of tran- 
sitions that contribute to the sum of Eq. |(SJ). Inter- 
branch transitions from the upper (A = 1) branch to the 
lower (A — —1) branch are possible everywhere inside the 
shadowed strips. They obey the oscillator selection rules 
(l,n) (l,n± 1). Intrabranch transitions are possible 
only in the immediate vicinity of the Fermi level, from 
the lower empty orbit to the upper filled orbit. 

The explicit expression for the contribution of intcr- 



branch transitions to Y^zxy is 



^zxy 



eh 



mai 



■E7 



(c„ + 1^/2) + vn{ 



.) 



(9) 

where the terms with v = ±1 describe the contributions 
from (l,n) (l,n± 1) transitions, respectively. For 
weak magnetic fields, when 7,rt ^ 1, the separate terms 
of the sum take asymptotic values l/2^/n, and summa- 
tion can be replaced by integration over n. Phase volume 
arguments suggest the relation n± — Pkj./2 between the 
integration boundaries n± and the Fermi momenta k± of 
the spectrum branches. Finally 

eh f "^ dn eh{k^ ^ ^'+) e 
I Snma Anh 



8TT\/2mal 



(10) 



note that two intraband terms T,^i:„ could cancel the sum 
of Eq. lO because they are large as n^/^ while the sepa- 
rate terms of the sum are small as n~-^/^. 

The above derivation of the zero spin-conductivity, 
'^zxy = 0, is valid in the uj/lOc 0,ujc limit and does 
not depend on any assumptions related to the scattering 
mechanism (provided the scattering is weak) and is ap- 
plicable for an arbitrary ratio of the spin-orbit coupling 
energy rac? jh} to [i. This derivation of the cancelation 
theorem provides a direct connection to the interbranch 
spin conductivity of Eq. H1U|) that has been found and 
generalized in a number of papers i^ii^^iSiiLiSiiSiSSiSS In 
what follows we provide a formal proof that shows that 
this theorem is valid at arbitrary B and follows from two 
sum rules. 

The first sum rule follows from the fact that the ma- 
trices {Xn\kx\X!ri!) and {Xn\ax\X!n') are symmetric and 
antisymmetric, respectively. Therefore, the numerator of 
Eq. ((SJ can be rewritten through a commutator of two 
operators, kx and Ox-, that commute. For a sum over A'n' 
(restricted only by the oscillator selection rule n' = 7i±l) 
we arrive at the identity 

^(An|L|A'n')(A'nV,|An) = i(An|[L, (T,]_|An) = 

A'n' 

(14) 
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that is valid for arbitrary {X,n). The second identity 
reads 

"^{Xnlk^lX'n + X){X'n + X\a.j,\Xn) = -2n/c„. (15) 



It is important that this sum does not depend on A. The 
latter identity can be checked by inspection using the 
eigenf unctions of Eq. Q . 

To calculate Y^zxy of Eq. ©, we introduce an auxiliary 
sum 



with n running across the shadowed strips of Fig. 1 and 
arbitrary A' and n' . The interbranch part of the sum, 
A' = —1, differs from 'S^.y only by two extra terms with 
n' near The intrabranch part of the sum, A' = 1, 
highly simplifies because the summand in Eq. IjlGI) is an- 
tisymmetric in the quantum numbers. In and In'. As 
a result, all terms but two, the last on the left and on 
the right, cancel. Finally, the equation relating ^zxy and 
^°xy reads 



^zxy — 



2-Kraa 



J2 J2^ln\kx\X'n'){X'n'\ax\lr 



n+ <n<n_ X'n' 



(16) 



J 



^zxy = Kxy + ^^[^^'^- ^ l|fcx|ln-)(ln-|(T,|ln_ - 1) + (ln_|fc,|ln_ + l)(ln_ + l|a,|ln_) 

+ (ln_|fc^|ln_ + l)(ln_ + l|a:,|ln_) - (In_|fc:,|In_ - l)(ln_ - l|CT^|In_) ) . (17) 



r 



In the parentheses, the two first terms arc interbranch 
contributions and the two last terms are intrabranch con- 
tributions. 

Now we are in a position to apply the sum rules of 
Eqs. lO and to Eq. (HTj). vanishes due to 

Eq. p4|l . The second and third terms in the parentheses 
can be summed by using Eq. 115|l . After transforming 
the first term by using the symmetry properties of the 
matrices kx and ax, the first and the fourth terms in 
the parentheses can be also summed, and the two sums 
cancel. Finally, we arrive at the cancelation theorem, 
'^zxy — 0, valid for arbitrary B. 

We conclude that the cancelation theorem of Eq. 
is a corollary of the sum rules of Eq. and It 

is an exact property of the free-electron Hamiltonian H, 
and spin conductivity should increase gradually with in- 
creasing electron scattering as has been found numeri- 
cally in Ref. E Therefore, dc spin conductivity of the 
systems described by the Hamiltonian of Eq. |Q is dissi- 
pative in agreement with the conclusions of Refs. .7.8L .9>.lCL 
The differences in the results of independent numerical 
g-g-Qj.^^io,23j24j25 g^j.g j^^^ properly understood yet and can 
originate from the differences in boundary conditions. 

It is instructive to follow up how the cancelation rule 
of Eq. H13() breaks down in ac regime. When lu ^ 0, each 
term in Eqs. lO and Hll|l should be multiplied by a factor 



trary, intrabranch terms E^^y(ti') change strongly at the 



XnX' 



(wAn — ^A'ri')^ 



(18) 



When u! <C 2ak±/h, these factors are close to unity 
for interbranch transitions, hence, the finite-frequency 
correction to "^Ixy can be disregarded. On the con- 



scale of LUc- When lu 
to the sum S^J.^^ 



zxy 



and 

quency region Uc 



^ ujc, the relative corrections 
are about SwV'^c- The sum 
diverges at w ~ ujc, changes sign. 



yll 

zxy 



zxy 

decreases as 



{uc/^)'^ for LU 3> ujc- In the fre- 
<C a; <C 2ak±/h, spin conductivity 
'^zxy{'^) ~ ^zxi/! hence, ac spin conductivity approaches 
the universal valued This conclusion correlates with the 
results by Inoue et alX and Mishchenko et al.^ who have 
found that in clean samples spin conductivity increases 
at cj ^ h/Tp, Tp being the momentum relaxation time. 
All these results emphasize the importance of transients 
for generating electrically driven spin currents 

The above results have been derived for a homogeneous 
system in a homogeneous electric field. Inhoniogeneities 
facilitate spin currents (and spin injection) both in the 
diffusiveSiSi and ballistic regimesii2i2& Mishchenko et al& 
have predicted developing transport spin currents near 
the electrodes, within a spin diffusion length from them, 
and spin accumulation near the corners. On the con- 
trary, equilibrium edge spin currents studied by Reynoso 
et alm^ are background currentsSS and do not produce 
spin accumulation. Spin currents near the free edges par- 
allel to the field £^ || y have not been investigated yet. In 
general, establishing a connection between spin currents 
and spin accumulation is a demanding and important 
problem. Indeed, despite the fact that transport spin 
currents can, in principle, be measured (because any de- 
viation of an electronic system from equilibrium produces 
electro-motive forces^) , and some procedures have been 
already proposedfSi it is the spin magnetization that is 
the principal measurable quantity. 
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The cancelation theorem of Eq. H13|l is vahd also for 
a Dresselhaus HamiltonianSi24 because it differs from 
the Rashba Hamiltonian only by an unitary transforma- 
tion (Jx <7y,ay <yx,(7z —Cz that does not change 
Eq. ©. 

For comparison, the Drude formula for Hall conductiv- 
ity (Jxy = {e'^niLjcTp/m)/{l + uj^t"^) has different ^ 

and limits, Cxy — and axy — e'^n/mujc, re- 

spectively. On the contrary, the diagrammatic {lOc = 0) 
and our (Tp^ = 0) results for dc spin-Hall conductivity 
match smoothly, what indicates a remarkable stability of 
its vanishing. 

In conclusion, we have shown that the cancelation of 



the inter- and intrabranch contributions to the dc bulk 
spin conductivity Jzxy is an intrinsic property of the free 
electron Hamiltonian and originates from two sum rules 
inherent in it. Weak electron scattering is only needed to 
reveal this cancelation in terms of the diagrammatic tech- 
nique. The cancelation persists in a magnetic field if the 
Zeeman interaction is disregarded, and spin conductivity 
increases resonantly at the cyclotron frequency. 
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Figure Captions 

FIG. 1. Schematic of the energy spectrum and elec- 
tronic transitions in a magnetic field. Here A = ±1 nu- 
merates spectrum branches, fc+ and k- are Fermi mo- 
menta for the upper and lower branches, respectively, n 
are Landau quantum numbers for these branches, n± are 
the lowest unoccupied states at both branches, and /j, 
is the Fermi level. Shadowed strips indicate regions of 
interbranch transitions between the occupied (— )-states 
and unoccupied (-l-)-states. The arrows n ^ n ± 1 show 
typical interbranch transitions in the "bulk" of the Fermi 
sea, while the arrows —t rt+ — 1 and ?i_ —t — 1 
indicate two intrabranch transitions at the Fermi edge. 
Cancelation of the inter- and intrabranch contributions 
results in the zero bulk spin conductivity Jzxy 




-3-2-10 1 2 3 



Momentum 



